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Abstract 

Quantum mechanical transition amplitudes are calculated within the stochastic quantiza- 
tion scheme for the free nonrelativistic particle, the harmonic oscillator and the nonrelativistic 
particle in a constant magnetic field; we close with free Grassmann quantum mechanics. 
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1 Introduction 

Stochastic quantization [1] was introduced more than a decade ago as a novel method for 
quantizing field theories and provides an intriguing connection between statistical mechanics 
and quantum field theory [2-3]. 

It is the aim of this paper to explore the somewhat more elementary issue of how to calculate 
within the stochastic scheme the transition amplitudes of ordinary quantum mechanics. Amus- 
ingly this problem has passed practically unnoticed over the years (see however [4-5]) so we 
would like to present it in a rather elementary way: We study the nonrelativistic free particle, 
the harmonic oscillator and the nonrelativistic particle in a constant magnetic field, presenting 
a phase space as well as a configuration space formulation; we close with the free Grassmann 
quantum mechanics [6]. 

At the beginning we recall basic facts on stochastic quantization: Given some arbitrary 
functional / of coordinates x(t), its vacuum expectation value can be obtained as the equilibrium 
limit s —* oo of the stochastic correlation function {f[x; s]) 

(f[x}) = l\m(f[x;s]) (1.1) 

which can be defined as 



(/[*;*]) = J Vxf[x]P[x;s}. (1.2) 



Here P[x; s] is a normalized (generalized) probability distribution, obeying the Fokker-Planck 
equation 

^-P[x;s} = I" * * ( ' (1.3) 



ds [X '' ai ~ Loo 6x(t) \5x(t) l 5x{t)) 

We work in real-time quantum mechanics and implicitly assume a convenient ie-insertion in 
the action S to guarantee the convergence of the stochastic process^ [7-10]. Alternatively the 
stochastic correlations (1.2) can be obtained by assigning to x(t) an additional stochastic time 
dependence x(t) — > x(t, s) via the Langevin equation 

dx(t,s) . 8S , , ,„ A . 

Its solution has to be inserted into the given functional / and the average over the white noise 
rj, characterized by 

{ri(t, s)rj(t', s')) = 26{t - t')5(s - s') (1.5) 

finally has to be worked out. 

In quantum mechanics in contrast to the calculation of vacuum expectation values one fre- 
quently is interested in the normalized transition elements of time ordered products of operators 
between Heisenberg state vectors \x,t) which are eigenstates of Heisenberg operators x(t) be- 
longing to the eigenvalues x. We therefore introduce new stochastic correlation functions, taking 
care of the boundary values Xi at ij and Xf at tt. We consider for example 

{xf,tf\T{x{ti)x{t 2 ))\x i ,t i ) 1(4 \ u ^ / r R \ 

, -r = hm {x(t 1 ,s)x(t2,s) ) bc (1.6) 

(Xf,tf\Xi,ti) s^oo 



1 Actually we would like to avoid a general discussion of complex stochastic processes but we note that this 
prescription is sufficient for the case of the quadratic systems under consideration. 
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and define in analogy to (1.2) 



(x(ti,s)x(t 2 ,s)) bc = Vx x(ti)x(t 2 )P[x;s]. (1.7) 

Jx(ti)=Xi,x(tf)=Xf 

P now satisfies a Fokker-Planck equation (1.3) with the time integration restricted to the 
interval t 6 [£$,£/]; the normalization condition reads (1)& C = 1. It follows that the stationary 
solution is modulo given by e lS so that the standard path integral representation emerges at 
the equilibrium limit of (1.7). 

Throughout this paper we prefer to use a Langevin equation approach which seems 
favourable given the rich technical experience already acquired in the case of field theories 
[2-3]; we will give the details of how to implement the boundary conditions in the Langevin 
scheme in the next section. Furthermore our ultimate goal (not yet addressed in this paper) 
will be the discussion of constrained quantum mechanical systems and we expect intriguing 
advantages of the Langevin approach over the rather complicated Fokker-Planck formulation 
[4]- 

The central object on which we would like to concentrate, however, is the transition ampli- 
tude (xf,tf\xi,ti) itself, which is one of the fundamental quantities in quantum theory. Due to 
the normalization condition for the probability density P we cannot directly express the transi- 
tion amplitude in terms of it. Stochastic expectation values calculated either with Fokker-Planck 
or Langevin equation techniques are always normalized automatically and it seemed difficult to 
find a way to reproduce such quantities as transition amplitudes within the framework of the 
stochastic quantization. It is, however, indeed possible to relate (xf,tf\xi,ti) to the normalized 
stochastic expectation value of the Hamiltonian H. This follows, as in the case of a regular 
quantum mechanical system the Heisenberg state vectors fulfill 

\xi,U) = Texpfi t H{t)dt}\ Xi ,t ) (1.8) 



t 



and 



Therefore we get 



d_ 

dti 



d 



ln(x f ,tf\xi,ti) = 



.(xf,tf\H(ti)\xi,ti) 



(1.9) 



(1.10) 



dU (xf,tf\xi,ti) 
The right hand side is a normalized expectation value of the Hamiltonian at U, which implies 
that this quantity is obtainable as a limit of the stochastic average (H(ti, s)) bc . Here note that for 
conservative systems the stochastic average of the Hamiltonian is i-independent at equilibrium 
s = oo; for calculational simplicity we evaluate it at the initial time t = ti. That is, we have as 
a first step 

rti 

(1.11) 



lim {H(ti,s)) bc dti\ 
s— »oo 



with c a constant independent of ti. For a transparent presentation we restrict ourselves to 
Hamiltonians quadratic in the canonical variables. In this case we can separate off very simply 
contributions from solutions to the classical equations of motion and arrive at 



i exp / / H cX {ti)dti 
cexp[i5 c i] exp 



exp 



\\m_{HQ{ti,s)) bc dti 
lim (H Q (ti,s)) bc dt. 



(1.12) 
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Here the classical action S c \ has appeared as a result of the famous relation dS c \/dti = H c \(ti) 
[11] and is composed of the classical path x c \(t) as usual. We remark that all the quantum 
contributions are contained in the second exponential factor. 

The constant c in principle could depend onxj and Xi which, however, is not the case. This 
follows from a similar variational principle as above for (xf,tf\xi,ti) w.r.t. Xf and x-i and the 
relations dS c \/dxf = p c \(tf) and dS c \/dxi = —p c \(ti)- So the constant c is indeed tf, ti, xj, 
Xi independent and can be fixed by requiring the transition amplitude to approach a Dirac 
5-function 5{xj — Xi) in the limit of ti = tf (compare also with Ref. [11]). 



2 Phase Space Formulation 

2.1 General Procedure 

In the last section we related the transition amplitude to an expectation value of the Hamilto- 
nian, so it appears appropriate to rely on a phase space formulation of stochastic quantization 
[12]. We put 

x(t) = x cl (t)+x Q (t), P(t) = Pd(t) + PQ (t) (2.1) 

and implement the boundary conditions for the quantum fields XQ^ti) = XQ(tf) = by the 
Fourier decompositions 

oo 

mr 



x Q (t) = ^ Xn sin— (t-U), (2.2) 

n=l 

oo 

pqV) = E^ cos Y"(*-^) + Y' T = t f -U- (2-3) 

71=1 

Notice that no boundary conditions are imposed on the momentum variable P Q(t) and that any 
function defined in [ti, U + T] can continuously be extended to [U — T,ti + T] as an even function 
like the above P q. Stochastic quantization proceeds by introducing the s-dependence for the 
Fourier modes x n — > x n (s),p n — > p n (s) according to the phase-space Langevin equations 

d .$Sq . d .5Sn , n 

— X n = l— + i n , —p n = + T]n (2.4) 

as ox n as op n 

where the noises fulfill 



(in(s)im(s')) = (Vn(s)Vm(s')) = 25 nm 5(s - s'), 

(Vn(s)U(s')} = 0. (2.5) 



For the quadratic case (2.4) is explicitly solved by 



x 



" )(,)=/ e iA ^ ) ( l n ) (a)da (2.6) 

y Pn J JO \ 'In ) 

where A n is a model-dependent matrix. Working in Stratonovich calculus we easily derive 



J nmi 



(z„(s)£m(s)) = (Pn(s)f] m (s)) = 5„ 
{X n (s)r]m(s)) = {Pn(s)€m(s)) = (2.7) 
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and find that correlations containing stochastic time derivatives vanish in the equilibrium limit 

s — ► oo 



x n (s) 



dx m (s) 
ds 



x n (s) 



dp m (s) 

ds 



. . dx m (s) 

Pn{s)- 



ds 



Pn(s) 



dpm(s) 

ds 



0. 



(2. 



In the following examples, we avoid solving the Langevin equations explicitly and instead extract 
the two-point correlations indirectly by using (2.8). 

2.2 Free Particle 

The easiest example to discuss is the nonrelativistic free particle with the Hamiltonian 



H 



P 

2M' 



We have now 



S, 



cl 



M 
2T 



(x f -Xi) 2 , S Q = J2 

n=l 



nir 



Pr, 



T 2 T 

T XnPn ~ 2M j 2 ~ Po 8M 



(2.9) 



(2.10) 



and get 



dx n 
ds 

dp n 
ds 

dpp 
ds 



2 Vn + snj 

/ nir p n \T 



' \ I U) 2 ' 



(2.11) 



From (p m (s) 
rium limit 



t d x n (s) 

ds 



Pm{s) - \ and (po{s) - ) we immediately find in the equilib- 



ds 



ds 



(PrnPn) = ^hn (p m (s)p n (s)) = 0, (PmPo) = 0, (pi) = 

respectively. Furthermore, 

{PQ{U,s)) bc _ (pi) i 



MM 



lim (Hn(ti,s))bc = lim 



s^oo 2M 



8M 2T' 



so that 



and finally 



■ |*\lhn (Hq^s))^ = -i J T (-^)dT = -hnT 



{x f ,t f \xi,ti) = c^=e lS -\ c=\\^- 



M 



2-Ki 



(2.12) 

(2.13) 
(2.14) 

(2.15) 
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2.3 Harmonic Oscillator 



Another standard example is given by the harmonic oscillator 

.2 1 



in which case 



Sq 



loM 



2 sin loT 

oo 



H = ^— + t-Mu?x 2 
2M 2 



[{xf + Xj) cos ljT — 2x{Xf] , 



so that 



71=1 

dx n 
ds 

dpn 

ds 
dPo 
ds 



Y. ( -f-''"/'" 77T7 ^ M ^l I - - Pi) 



1 



2M 2 



"/ 2 



/ 2 _T_ 

8M 



n7r 



i I — 1„ - — ) — + 77„, 



It'" m; 2 

- i lM P0 + r >°- 



From {p n ~j~^J i W ™ ^ e ec l um ^ rmm limit 

. . . 2 . 4M 

(PnPo) = 0, (po) = -t-=- 



and from ( p m —r^-\ combined with ip m -^— 
ds \ ds 



{PnPm) = 



2iM 



.Too) 



Therefore (recalling the boundary condition xq(^) = 0) 

hm (HQiU,*))* = ( E (PnPm) + ^ 



2M . 

\ n,m=l 

and the transition amplitude is obtained as 

1 



— cot 
2 



JS C ] 



V 7 sin Tcj 



c = 



2tt7' 



(2.16) 

(2.17) 
(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



2.4 Constant Magnetic Field 

A less trivial example of our approach is the nonrelativistic charged particle in a constant 
magnetic field: 



H= l -{p-Af, 



P 



( p \ 




( o > 






q 


Bx 


- 




V r J 




\ o } 







(2.24) 



We split again 



and obtain 



x(t) = x c i(i) +x Q (t), p(t) = p c \(t) + p Q (t) 



1 \(z f -Zi) 2 B BT,, x2 , ,21 „/ 



as well as 



T _ 



n=l 



£ 5 



£T 

7T 

The Langevin equations read as 



(-l) n -l ^ (-i)n+m_ 1 

m=l 



2n 



nir 



~ds~ = 

dp n . (_ rnr _ T j> 

Po 
ds 



1 ( - jPo + c) +r] 



where we have introduced 



BT 



a n = 








\ 



K = I /?n 







< - ; 7 




with 



a n = 



7 



(-1)™ - 1 2 ~ _ j 

-5x n - <? h - > g m n = — 

nir ir ^-r. n z — m z 

m=l 



^ (-l)n-Hn _ ! 
Pre — / , re 9 o 

m=l 



ST 



2vr ^ 



(_l)m _ j 



m=l 



It follows from ( qvn^r 1 ), ( go-T^Y (pn-zr-\ ( r m^T L \ and (r n ^- ) that 



(is 



7 / 7 \ Jr it> 7 / ? \ tit 7 / -~ \ • 1 1 7 

as I \ as I \ as I \ as 
(q n q m ) = (q n qo) = (PnPo) = (r n r m ) = (r„r ) = 



Therefore 



1 1 00 

Hm {H Q (U,s)) bc = -(po 2 ) + - (PnPm)- 



n,m=l 
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Using (p®~lf^^> ^\ r °~d~^ S } W6 ^ ave as P rev i° us ly 



From / qo ■ 



I dq \ 
\ q °~ds) WG 



(2.32) 

we relate (qop n ) to (x n qo); then (^o - ^ 1 ^ connects (x n qo) to ((j>g) and from 



From (pm-r 11 ) we 



we obtain 



<gg) = -i25cot— . 
-j^j we relate {x n p m ) to {p n Pm); from (prn^^j relating {x n p m ), {q p n ) and 



(2.33) 

,2\ 



B 2 T 

(PnPn>) = 2« 2~ 



L 2 ~ 

Straightforwardly we finally arrive at 



S nn , BT (-!)"-! (-!)"'-! BT 



and 



, itt i B BT i 

hm {Hq(U, s)) bc = -i— cot — - — 



(Xf,tf\Xi,U) 



= c 



2 UvriJ 



3/2 



Tsin^ 



(2.34) 

(2.35) 
(2.36) 



3 Configuration Space Formulation 

3.1 General Procedure 

Although as outlined in the last section the phase space formulation is the natural one, it is for 
technical reasons preferable to rely on the simpler configuration space formulation. If 

~2 M _ 

(3.1) 



H = ^- + V(x), C = —x 2 -V{x) 
2M y h 2 w 



we have, of course, 



(JI(?(I ,. s)) hl - (p Q (t u s)^jQ(ti,s)^ - {£ Q (U,s)) bc . 



Now 



{C Q (ti,s)) bc 



M ldx 2 



dt 



M 



be 



— lim d tl d t2 {xQ(t 1 ,s)xQ(h,s)) bc 

Z tl,t2— *H 



(3.2) 



(3.3) 



where we have used Feynman's time splitting procedure to properly define the kinetic energy 
contributions [11]. The transition to a pure configuration space formulation is achieved by 
observing that generally 



dpQ 
ds 



(t,s) = i 



( dx Q (t,s) _ OHq \ 

{ dt d PQ (t,s)) + ^ s) - 



(3.4) 
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From / — — — (t' s) (t,s) ) we find 
at os 



be 



lim 

s^oo 



(3.5) 



and finally obtain 



M 

lim (H Q (ti,s)) bc = — lim d tl d t2 lim {x Q (t u s)x Q {t 2 , s)) bc , (3.6) 

S — >00 ^ Cl,*2 — >H s — 

which coincides with the very definition of the Hamiltonian in the path-integral formulation 
[11] and can be calculated entirely in configuration space. We start from the regular Lagrangian 
C and separate the coordinate x(t) into a classical and quantum part, the latter being Fourier 
expanded as in (2.2). The configuration-space Langevin equations for the Fourier modes are 
associated as in the original Parisi-Wu approach and the transition amplitudes can be calculated 
using (1.12) and (3.6). We remark that the transition to configuration space in the case of 
the particle in a constant magnetic field arrives at the same expression (3.6), even though it 
corresponds to a case of velocity-dependent quantum mechanical potential. 

3.2 Free Particle 

To demonstrate the configuration space approach we start from the free particle action 



S = J t dt—x , S Q = }_^ 4T x n , (3.7) 

and get 



so that 



From this we immediately find 



dx n .nVM 
2iT 

\X n X m ) = t, 9 &nm ■ (3-9) 



1 

}^n o {xQ(t 1 ,s)x Q (t 2 ,s)) bc = j^;[tf - max(ti,t 2 )][min(ti,t 2 ) ~ U] (3.10) 
and have from (3.6) 

^ o {H Q (U,s))bc = ~- (3-11) 

Here we have to keep in mind that we should first sum over Fourier modes n to obtain 
lim^oo (xq (ti , s)xQ(t2, s)) bc = {xQ{t\)xQ(t2))bc-> differentiate w.r.t. t\ and t 2 and finally take 
the limit ti,t 2 — > t%. This ordering corresponds to Feynman's time splitting procedure [11], 
which avoids infinite contributions arising from the velocity correlation at the same point. If 
the order were reversed or changed, for example, if we differentiated xq{€) before the summation 
over n, we would have an infinite contribution Y^=i = °°; as i s easily seen from (3.9). 
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3.3 Harmonic Oscillator 

With 



S = jf ' dt (fx 2 - , S Q = £ 



4T 



71 



V 7T 



x n> 



we obtain now 



and 



dx„ . Mir 2 



ds 2T 



2 f T ^\ 
n - I — 1 x„ + r/ n 



/ \- ■ 2r 
\x n x m ; - « M7r2 



Furthermore we obtain 
and get upon differentiation 



lim (Fq(^, s)) 6c = — -cot Tuj. 



3.4 Constant Magnetic Field 

Despite the simple looking form 

Sq takes an involved expression 





oo 


/ 


D 


M 






' x n > \ 






-M 


D 






Vn> 


n,n'=l 












nn' 





where 



Dnn' — fl> $nn' i 
M n „f = 



(-l) n+n '-l ,2BT 
nn 



n 2 — n' 2 



7T 



2 ' 



(3.12) 

(3.13) 
(3.14) 



smuj(tf — max(ti, t2)) sina;(mm(ii, £2) — U) (3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



As the ^-components decouple, we first investigate the (x n ,y n ) contributions only. We start 
from the Langevin equations 



d_ ( x n \ _ijf_ ( D M 
~ds \ Un j ~ 2T 1 —M D 



Cn 



(3.20) 



and multiply with (%» , y n n ) to extract in the equilibrium limit 



I D M 
2T [ -M D 



(x n ",y n ") 



Inn" 



(3.21) 
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We eliminate the off-diagonal expectation values and find for x n (we remark that an identical 
expression holds for y n ) 



- % -^{D + MD 1 M) nnl {x n ,x n „) = <W' 
which explicitly evaluates to 

.2 rr /d^\21 /BT\ 2 [(-1)"- i][(-i)"' - 1]] 



ITT 

2T 



n 



(f) 2 



S nn > + 2 



V vr 2 J 



nn' 



(3.22) 



(x n 'X n ") = $nn"- (3.23) 



We find 



\X n X n r) g 



(BTf [(_!)»_ !][(_!)»' _i] 



7T^ 



n- 



2 Tr" 



nn' 



n 1 - 



cot 



(3.24) 



Therefore after the summations over n and n', we get the correlation function 



(xQ(ti)x Q (t 2 ))bc = 



;{sinB(tf — max(ti,t2))sinS(min(ti,t2) — t-i) 



B sin BT 1 

B BT B BT 

-[cos— (tf +U- 2ti) - cos— J[cos— (t f + U - 2t 2 ) - cos— J), 



from which follows 



iB BT 

lim d h d t2 lim (x Q (ti, s)x Q (t 2 , s)) bc = — —cot—. 



(3.25) 



(3.26) 



Adding an identical ^-contribution and the trivial ^-contribution we get in total as before 



.. iB BT i 

hm (H(ti, s))bc = cot — . 

s-»oo x v u " oc 2 2 2T 



(3.27) 



3.5 Nonrelativistic Grassmann Quantum Mechanics 



In this section we want to extend our method to deal with Grassmann quantum systems. We 
study the simplest possible model, which is a free nonrelativistic Grassmann particle described 
by [6] 



£> = — , 

2 V dt J 

Here we have introduced a Grassmann vector 

= {e a }, a = 1,2 

and defined an inner product 



e 2 = e • e = e Q r Q/3 e /3 



with the help of a metric tensor 



r a/3 = -r°0 



o 1 
-1 o 



(3.28) 

(3.29) 
(3.30) 

(3.31) 
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Notice that (3.28) differs from Grassmann quantum mechanical models which are inspired from 
particle physics and have kinetic terms linear in the time derivative [13-15]. These models 
generally lead to second class constraints, however. Having decided to restrict ourselves in this 
paper to regular systems we do not include these models in our discussion. 

We choose to work in the configuration space formulation and separate the Grassmann 
variable into a classical and quantum part as before 



mr 

W n sin 

n=l 

so that 



e(t) = e cl (t) + e Q (t), e Q (t) = £e„sin— (t-u) (3.32) 

=1 

Sol = fl5^M S Q -±^i. ,3,3) 

n=l 

The definition of the Langevin equation for the Grassmann mode 0„ requires (similarly to the 
fermionic field theory case [16]) the introduction of a kernel, which we choose to be F al3 

d@ n _ ;-na0 SS Q ■ _q 

= ^<S + C (3.34) 

KMis')) = 2S nm T^S(s - s') (3.35) 

ds 



Here the noise fulfills 

which guarantees the Grassmann anticommutativity. From (@ m (s) ■ ^f" (s) ) we find 



AiT 

lim (@ n (s) ■ Q m (s)) = , o.J nm (3.36) 

where a crucial minus sign has appeared as a consequence of the definition of the inner product 
with T a/3 . We easily deduce 

2i 

}}*% ( Q Q( tl > s ) " ®Q(*2,s))b c = -Jjfltf ~ max(ti,t 2 )][min(ti,t 2 ) - U] (3.37) 
so that after differentiation 

]im o {H Q {t i ,8))bc = ^ (3-38) 
and the desired transition amplitude becomes 

(e f t f \e i t i )=cTe iS «, c=^. (3.39) 

Note that the constant c was fixed in this example by requiring the transition amplitude to 
approach a Berezin delta function in the equal time limit t% = tf [6] . 
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4 Conclusions 

In this paper we have filled a somewhat historical gap in the applications of stochastic quan- 
tization and calculated for the nonrelativistic free particle, the harmonic oscillator, the non- 
relativistic particle in a constant magnetic field and the free Grassmann particle the quantum 
mechanical transition amplitudes. Our main observation is the possibility to relate the non- 
normalized transition amplitudes with the normalized exectation values within the stochastic 
quantization scheme taking care of the boundary conditions. 

Our procedure can equally well be applied to Euclidean quantum mechanics, in which case 
the corresponding partition functions become straightforwardly calculable. 

It remains a challenge to apply our scheme to constrained systems. 
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